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Abstract 

By extending previous analysis of the authors, a systematic study of the singularity 
structure and possible asymptotic behaviors of five-dimensional braneworld solu- 
tions is performed in the case where the bulk is a mixture of an analog of perfect 
fluid (with a density and pressure depending on the extra coordinate) and a mass- 
less scalar field. The two bulk components interact by exchanging energy so that 
the total energy is conserved. In a particular range of the interaction parameters, 
we find flat brane general solutions avoiding the singularity at finite distance from 
the brane, in the same region of the equation of state constant parameter 7 = P/ p 
that we found previously in the absence of the bulk scalar field (—1 < 7 < —1/2). 



^On leave from CPHT (UMR CNRS 7644) Ecole Poly technique, F-91128 Palaiseau, France. 



1 Introduction 



In our previous works [H |2], we started a systematic study of the singularity structure 
and possible asymptotic behaviors of five-dimensional braneworld solutions, by param- 
eterizing the bulk field content with an analog of perfect fluid satisfying the equation 
of state P = 7p, where the 'pressure' P and the 'density' p depend only on the extra 
dimension Y and 7 is a constant parameter. Our motivation was based on the idea of 
the so-called self-tuning mechanism for the cosmological constant ^ i4j, aiming to ex- 
amine in a model independent way the possibility of avoiding singularities in the bulk 
at a finite distance from the brane position. We had found three regions of 7 leading to 
qualitatively different behavior: 

• The region 7 > —1/2 is very similar to the case of a massless bulk scalar field. In- 
deed, the existence of a singularity at a finite distance is unavoidable in all solutions 
with a fiat brane, in agreement with earlier works that made similar investigations 
in different models, using other methods [5l [6]. Moreover, we have shown that 
the singularity can be avoided {e.g. moved at infinite distance) when the brane 
becomes curved, either positively or negatively. Thus, requiring absence of singu- 
larity brings back the cosmological constant problem, since the brane curvature 
depends on its tension that receives quartic divergent quantum corrections. 

• In the region — 1 < 7 < —1/2, the curved brane solution becomes singular while 
the fiat brane is regular. Thus, this region seems to avoid the main obstruction of 
the self-tuning proposal: any value of the brane tension is absorbed in the solution 
and the brane remains fiatH 

• In the region 7 < —1, corresponding to the analog of a phantom equation of state, 

the brane can be ripped apart in as much the same way as in a big rip singularity. 

^Of course, the main question is then whether there is a consistent field theory reahzation of such a 
fluid producing naturafly an efl'ective equation of state of this type [7] ■ 
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This happens only in the flat case, while curved brane solutions develop 'standard' 
collapse singularities. No regular solution was found in this region. 

Moreover, we have shown that all possible singularities at a finite distance from the 
original position of the brane can be classified in three main classes which we coined 
collapse type I, collapse type II, and big rip singularities, respectively. 

The collapse types I and II met in the asymptotic evolution are both characterized, 
as their name suggests, by the vanishing of the warp factor. Their differences can be 
traced in the behavior of the derivative of the warp factor and density of the matter 
component. In the collapse type I class for instance, the derivative of the warp factor 
becomes infinite whereas in the collapse type II class it remains finite. The density of the 
bulk matter, on the other hand, is necessarily divergent asymptotically in the collapse 
type I class, whereas in the collapse type II class it may approach a constant or even 
vanish. 

An interesting aspect of these two behaviors is that the types of singularity which 
become asymptotically feasible depend in the first place, on the type of bulk matter: 
while a massless scalar field, which may be regarded as a fluid with 7=1, allows the 
development of only a collapse type I singularity, a perfect fluid allows in addition the 
emergence of singularities covering the whole variety of the collapse type II class as well 
as big rip singularities. The latter are singularities characterized by the divergence of the 
warp factor, its derivative and the matter density of the fluid, and arise only when the 
parameter 7 is less than —1. In addition, the possible types are determined by the spatial 
geometry of the brane: A flat brane allows the development of all the different types of 
finite-distance singularity whereas a curved brane permits exclusively the formation of 
collapse type II singularities. 

In this paper we extend our previous work by examining the case of a bulk filled 
with a mixture of fluid and massless scalar field. We let these two bulk entities either 
interact with each other or simply coexist independently in the bulk. In the latter case, 
we show that all previous types of singularity are still possible. However, the nature of 
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each type is now enriched with the behavior of both bulk components. The collapse type 
I singularity, for example, is characterized by the divergence of the density of either one 
or even both of the two components, whereas, the big rip singularity is characterized 
by the divergence only of the density of the fluid while the density of the scalar fleld 
vanishes asymptotically. These two types of singularity arise in general solutions. The 
collapse type II singularity on the other hand, arises in particular solutions and exhibits 
a possible divergence in the density of the fluid. Apart from these singular solutions, we 
also flnd regular ones that lead to avoidance of flnite-distance singularities but only for 
the case of curved branes. In particular, we may avoid flnite-distance singularities for a 
curved brane when 7 > — 1/2. In contrast with the previous results, we do not flnd any 
range of 7 that leads to avoidance of flnite-distance singularities for a flat brane. As we 
mentioned before, the existence of a general regular solution for a flat brane implies that 
a self-tuning mechanism may be constructed. The failure of our flat-brane model to offer 
such possibility leads us to its generalization which is implemented by considering an 
interaction between the two components in the bulk. In this more complicated case we 
flnd that for an adequate choice of the interaction parameters and for —1 < 7 < —1/2, 
the avoidance of singularities is recovered. 

Our approach in avoiding flnite-distance singularities is to flnd ranges of the param- 
eter 7 and later on of the interaction coupling-coefficients that allow for the existence of 
solutions that are singular only at infinite distance. We do not, however, consider the 
possibility of constructing regular solutions with a matching mechanism as in [7]. This 
mechanism is implemented by exploiting solutions that exhibit a finite-distance singu- 
larity that is located only in the half line of the extra dimension away from the position 
of the space of matching. 

The structure of this paper is the following: In Section 2, we start by giving a set up 
of the basic equations of our model consisting of a brane in a bulk with a scalar field and 
an analog of perfect fiuid. These two bulk components may exchange energy in a way 
that the total energy is conserved. The field equations are written as a dynamical system 
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which we analyze with the method of asymptotic sphttings, cf. [8]. As a first step in 
Section 3, we focus on the case in which there is no exchange of energy between the bulk 
components and derive all possible asymptotic decompositions of the dynamical system 
together with their dominant balances, i.e., the different possible asymptotic modes of 
behavior. In particular Subsections 3.4-3.6, are devoted to the asymptotic structure of 
our braneworlds near finite-distance singularities in the bulk, while in Subsection 3.7, 
we focus on behavior at infinity. As a second step in Section 4 we consider that the two 
bulk components interact with each other and resolve the unwanted situation discussed 
in Subsection 3.7. We conclude and discuss our results in Section 5. 



2 Field equations 

We consider a braneworld model consisting of a three-brane embedded in a five- dimensional 
bulk space that is filled with a massless scalar field and an analog of perfect fluid. We 
assume a bulk metric of the form 

g, = a\Y)g, + dY\ (2.1) 
where is the four-dimensional flat, de Sitter or anti de Sitter metric, i.e., 

g, = -dt" + /2^3, (2.2) 

where 

g^ = dr"" + hlg2 (2.3) 

and 

^2 = rfe2 + sin2(^)V- (2.4) 

Here = 1, cosh.{Ht) / H, cos{Ht) / H {H^^ is the de Sitter curvature radius) and = 
r, sinr, sinhr, respectively. For the scalar fleld we assume an energy-momentum tensor 
of the form = (pi + Pi)uaUb — PigAB where A, _B = 1, 2, 3, 4, 5, ua = (0, 0, 0, 0, 1) 
and pi. Pi is the density and pressure of the scalar fleld which we take as Pi = pi = 



4 



A0'^/2, where the prime denotes differentiation with respect to Y and A is a parameter. 
Respectively, the energy-momentum tensor of the fluid is = (p2 + P2)uaUb — P2gAB 
and we assume an equation of state of the form P2 = 7P2 between the pressure P2 and 
the density p2 with 7 being a parameter. All quantities pi, p2 and Pi, P2 are functions 
of the fifth dimension Y only. The five- dimensional Einstein field equations, 

Gab = 4Tab, (2.5) 
where = M^"^ and M5 is the five dimensional Planck mass, can then be written as 

^ = _AA0'2_^A(1 + 27)P2, (2.6) 
a 6 

_ = + + _ (2.7) 

where A = Kg/4, /c = 0, ±1 (and the prime (') denotes differentiation with respect to 
Y). We assume that there is an exchange of energy between the two matter components 
depending on the values and signs of the two constants u, a, such that the total energy 
is conserved so that we have the following two equations, 

, , ; , /; . , Ct' , /2 AZ^ q! , ,0 0! , , 

A0y + 4A-C = -- — C + ap2-, (2.8) 
a la a 

p'2 + 4(7 + l)-p2 = — -0''-ap2-. (2.9) 
a Z a a 

Eqs. ( 12. 6 p and (12. 7p are not independent, since Eq. (12. 6p was derived after substitution 

of Eq. (12. 7p in the field equation Gaa = i^^T^a = 4:ATaa, a = 1,2, 3, 4: 



a" a'^ kH^ 2A,^,2 4A 

- + — - ^ = -^M' - ^1P2. (2.10) 



In our analysis we use the independent Eqs. (12.61) . (12. 8p and (12. 9 p to determine the 
unknown variables a, a! , 0' and p2, while Eq. (12. 7p will play the role of a constraint 
equation for our system. 

Our purpose is to find all possible asymptotic behaviors of (general or particular) 
solutions of the system defined by the dynamical equations (I2.6p -( l2l9l) . The most ade- 
quate tool for this quest is perhaps the method of asymptotic splittings summarized in 
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[8] . The first step is to write this system in the form of a suitable dynamical system. We 
introduce the following set of variables: 

{x,y,z,w) = (a,a',0',p2). 

The system of equations (12. 6p . ( 12. 8 p and (12. 9p then becomes the following dynamical 
system 

x' = y (2.11) 

y' = -AXz^x--A{l + 2-f)wx (2.12) 

= -(4+^)fc + f^ (2.13) 

V 2/ X X xz 

u.' = -(4(,+ 1)+.)!^!^ + 4^H:. (2.14) 

X 2 X 



while equation ( 12. 7p now reads 



y^ AX 2 2A kH^ 



z' + —w + (2.15) 



Since this last equation does not contain derivatives with respect to F, it is a constraint 
equation for the system (I2.1ip - (l2.14p . The vector field defined by the above system is 
given by 

f ( /i\ 2 2 f u\yz ayw yw Xuyz^Y 

f = ^y, -AXz X - -Ail + 2^)wx, " (4 + - j - + --, -(4(7 + 1) + ^) ^ + T 

(2.16) 

Before we proceed with the analysis of the above system, we introduce the following 
terminology for the possible singularities to occur at a finite-distance from the brane. 
Specifically we call a state where: 

i) a — J- 0, a' — )■ 00, (j)' — )■ 00, P2 — ^ 0, ps, 00: a singularity of collapse type I, 

ii) a — )■ 0, a' — !■ a^, 0' — )■ 0, p2 — ^ Ps, 00: a singularity of collapse type II, 

iii) a— >oo, a'— 7-— 00, 0'— )-0, P2— ^00: a big rip singularity, 
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where a'^ and ps are non-vanishing constants. 

In the following Subsections we first analyze the case in which there is no exchange 
of energy between the two components in the bulk, that is we take u = a = 0, and later 
on we examine the very interesting case cr = and u arbitrary and we comment on the 
results of the case z/ = and a arbitrary. The generic case a, v nonzero is difficult to 
study in generality and classify all possible behaviors. 



3 Non-interacting mixture in the bulk 

In this Section we let = cr = so that the system Eqs. (I2.1ip - (l2.14p becomes 

x' = y (3.1) 
2 

y = -AXz^x - -A{1 + 2-f)wx (3.2) 

z' = -4^ (3.3) 

X 

w' = -4(7 + 1)^, (3.4) 

X 

while equation fl2.15p remains the same. The vector field of the above system is 



T 



f = h/,-AA2'x- -A(l + 27)w;x,-4^,-4(7 + l)^ , (3.5) 
\ 3 X X J 

and there are three possible ways of decomposing it. We analyze them in turn in the 
following Subsections. 

3.1 Decomposition I 

The first way of decomposing the vector field (13. 5p is to assume that its dominant part 
is given by the form 

^ T 

X ' ' ' ' X ' 

while its candidate subdominant part is: 



f(°) = (y, -AXz\ -4^, -4(7 + 1)^)^ , (3.6) 
V X X / 



f« = (^0,-^A(l + 27)«;x,0,0^ . (3.7) 



7 



At this point we wish to determine all dominant balances, that is pairs of the form 

S = {a,p}, where a = (a, /3, c, C), p^ip,q,r,s), (3.8) 

with 

(p,?,r,s)eQ^ and (a, /3, c, () G x {0}, (3.9) 

that describe all possible asymptotic behaviors around the assumed position of the sin- 
gularity at Yg. We thus insert 

(x, y, z, w) = (aTP, ^T^ cT^ C^, (3.10) 

where T — Y — Yg, into the asymptotic system defined by the first decomposition, that 
is 

(3.11) 
(3.12) 
(3.13) 
(3.14) 

For each balance we need to 

(3.15) 



X ^ y 

y' = -AXz^x 



X 

yw 



w' = -4(7 + 1) 



X 



This leads us to the list of all possible dominant balances, 
check that the dominance condition, 

f(i)(aTP) 
lim ^ , ' = 0, 
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is satisfied and then discard those balances that do not satisfy Eq. fl3.15p . We end up 
with the following acceptable balancej^: 

= {(«,«/4,y3/(4v^), 0,(1/4, -3/4,-1, -(7 + 1))}, 7 < 1, (3.16) 

jB2 = {(a,a/4,y3/(4v^),0),(l/4,-3/4,-l,s)}, (3.17) 

1B3 = {(a, a, 0,0, (1,0, -1,-4(7+1))}, 7 < -1/2, (3.18) 

jB^ = {(«,«, 0,0),(l,0,-l,s)}, (3.19) 

jB^ = {(a,0,0,0),(0,-l,-l,s)}. (3.20) 

The balance 1B5 leads to p2 being identically zero which means that it describes a be- 
havior that applies in the case of bulk filled exclusively with the scalar field. We studied 
this balance in our previous work in [2] and found it to be unacceptable. We therefore 
do not consider it any further in this paper. 

The balances /i3i_4 are exact solutions of the system (13.1ip - fl3.14p . The constraint 
Eq. f l2.15p has contributed in the Eq. f l3.12p with all its terms excluding the term of 
the fluid density. We can therefore substitute these balances in the constraint equation 
f l2.15p neglecting the term of the fluid density and find out if they correspond to a fiat or 
curved brane. We find that for 7 7^ —1/2, the balances jBi and 1B2 correspond to a fiat 
brane while the balances 1B3 and 1B4 correspond to a curved brane with the arbitrary 
constant a satisfying = kH^. 

The value 7 = —1/2 is of special interest for our analysis since for this value of 7 
the system (I3.ip - (13.4p identifies with the system of this first decomposition (I3.1ip -( l3.14p . 
This means that the balances we find for 7 = —1/2 are exact solutions of the system 
f l3.ip -f l37i|) and they should therefore satisfy the entire constraint equation fl2.15p . These 
balances are: 1B2 (flat brane), 1B4 (curved brane with = kH^), fB^ (fiat/curved brane 
with C = 3/(2A)(l - kH^a^) and jBi (curved brane with C = -3/{2A){kH^/a'^)). 



^In the balance iBi the coefficient V3/(4-\/A\) may also be —VS/ (4%/ AX). This is true for every 
balance we find that has a square root in the coefficient of (j)' , but for simplicity we examine these 
balances only for the (+) sign. 
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3.2 Decomposition II 

The second way of decomposing the vector field (13. 5p is to take its dominant part to be 

f^°)= fz/,-|A(l + 27)^x, -4^,-4(7 + (3.21) 
\ 3 X X J 

Its candidate subdominant part reads, 

i^'^ = iO,-A\z^x,0,Oy. (3.22) 
For this second decomposition the system is given by the following equations 

x' = y (3.23) 

y' = --A{l + 2-f)wx (3.24) 
3 

z' = -4^ (3.25) 

X 

w' = -4(7 + 1)^, (3.26) 

X 

and the acceptable balances are calculated to bqf 

nBi = {(a,ap,c,3pV(2A)),(p,p-l,-4p,-2)}, I7I > 1, (3.27) 

11B2 = {{a,ap,0,3py{2A)),{p,p-l,r,-2)}, 7^^-1,-1/2, (3.28) 

ijB, = {(«,«,0,0),(l,0,r,-2)}, 7 7^-1/2, (3.29) 

iiB, = {(a,0,0,0),(0,-l,r,-2)}, 7 7^-1/2, (3.30) 

where p = 1/(2(7 + !))• Following the same trend as we did for Decomposition I, 
we see that the constraint Eq. fl2.15p has contributed in the Eq. fl3.24p with all its 
terms excluding the term of the derivative of the scalar field. We therefore substitute 
the balances uBis in the constraint equation Eq. f l2.15p neglecting the term of the 
derivative of the scalar field and find that the balances uBi and 11B2 correspond to a 
fiat brane while the balance uB^ corresponds to a curved brane with = kH^. 



■^The balance 11B4 is not analyzed any further since it leads to 0' being identically zero and a similar 
argument applies as in the case of 1B5 discussed above, in the decomposition I. 
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3.3 Decomposition III 

The third way of decomposing the vector field (I3.5P is to assume that all terms are 
dominant so that the system is given by Eqs. fl3.1l) - fl3.4p . For this third decomposition 
the dominant balances are 

HiBi = {(a,a/4,c,3/(32A)-AcV2),(l/4,-3/4,-l,-2)}, 7 = 1, (3.31) 

jnB2 = {(a, ap,0,3pV(2A)),(p,p- 1,-1,-2)}, 7^-1,-1/2, (3.32) 

niBs = {(a,a/4,v^/(4v^),0),(l/4,-3/4,-l,-2)}, (3.33) 

HiBi = {(a, a, 0,0), (1,0, -1,-2)}, (3.34) 

with p = 1/(2(7 + !))• These balances are exact solutions of the system fl3.ip -f l37il) and 
they should therefore satisfy the constraint equation ( ]2.15p . We find that the balances 
iiiBi, 111B2 and niBj, correspond to a fiat brane, while the balance niBi corresponds to 
a curved brane with = kH^. Notice that the three decompositions considered above 
exhaust all possible asymptotic ways that the vector field (13. 5p can split. 

Subsections 3.4-3.7 are the heart of this Section that focuses on non- interacting bulk 
components. We have grouped the possible balances i-mB found above into four dif- 
ferent sets according to the type of singularity they lead to, or, their regular behavior. 
For each particular balance we follow the method of asymptotic splittings up to the 
point where we end up with a well-defined series expansion. These expansions serve to 
completely justify our claims that the asymptotic behavior of the braneworld is the one 
claimed. We find that in all cases these behaviors result from Puiseaux representations, 
in particular there are no logarithmic terms present in any of the expansions. 

3.4 Collapse type I singularities 

In this Section, we analyze the balances jBi, niBs, 1B2, 11B2, iiiBi and 111B2 that 
describe the asymptotics around collapse type I singularities. 
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3.4.1 The balance iBi 

We start with the analysis of the balance jBi that corresponds to a flat brane for 7 7^ 
— 1/2. We will show that for different values of 7 this balance implies different behaviors 
of the matter density of the fluid around a collapse type I singularity. We flrst have to 
calculate for this balance the /C-matrix given by 

//Ci = (a) - diag p, (3.35) 

where Df^^^a) is the Jacobian matrix of the dominant part f*^°^ in Eq. (13. 6p . 

/ 1 \ 



x,y,z,w) 



-\Az^ 



-2\Azx 



-4- 



_4i 



X 



X 



4(7 + 1)^ -4(1 + 7)- -4(1+7)^ ) 

X'' X X ' 



(3.36) 



evaluated on a. We have that a = (a, a/4, a/3/(4v^AX), Q and p = (1/4, —3/4, —1, —(7+ 
1)), so that the /C-matrix in this case is 

,/Ci = DfW((a, a/4, v^/(4yiA),C))-diag(l/4, -3/4,-1, -(7 + 1)) = 



1 

~4 
3 

~16 







3 

4 

%/3 



a 



4aA/AX aVAX 







I (1 + 7)- -4(1 + 7)- 

\ a 



a 
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Next, we calculate the /C-exponents for this balance. These exponents are the eigenvalues 
of the matrix //Ci and constitute its spectrum, spec(//Ci). We wish to build series 
expansions of the variables in the form 

X = TP(a + E°liCjT^/^), (3.37) 

where x = (x, z, w), Cj = (cji, Cj2, Cjs, Cji), and S is the least common multiple of the 
denominators of the positive /C-exponents and the non-dominant exponents q^^^ defined 
by the requirement 

, (3.38) 

(cf. 0, [in])- The arbitrary constants of any particular or general solution first appear 
in those terms in the series (13.371) whose coefficients have indices k = gS, where g is 
a non- negative /C-exponent. The number of non- negative /C-exponents therefore equals 
the number of arbitrary constants that appear in the series expansions of (I3.37p . There 
is always the —1 exponent that corresponds to an arbitrary constant that is the position 
of the singularity, Yg. 

The balance jBi corresponds thus to a general solution in our case if and only if it 
possesses three non-negative /C-exponents (the fourth arbitrary constant is the position 
of the singularity, Yg). Actually, once we use the constraint equation (I2.15p . one of the 
three arbitrary constants corresponding to the three non-negative /C-exponents will be 
set to a specific value, so that the general solution of the dynamical system ( ]3.ip -( !334l) 
with constraint equation (12.151) will exhibit three in total arbitrary constants taken into 
account also the singularity position Yg. Here we find 

spec(//Ci) = {-1, 0, 0, 3/2}. (3.39) 

The double multiplicity of the zero /C-exponent refiects the fact that there are two ar- 
bitrary constants in this dominant balance. In total we have three non-negative /C- 
exponents which means that this balance indeed corresponds to a general solution. 

Naturally, the behavior of p2 depends on the exponent —(7 + 1). We try inserting 
different values of 7 so that to trace all possible asymptotics for p2 keeping in mind that 
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7 < 1. For instance, for 7 = we have that 

iB, = {(«, a/4, V3/{4VAX), (), (1/4, -3/4, -1, -1)} (3.40) 
and substituting in the system fl3.ip -( l3^ the particular value 7 = and the forms 

— yoo y-j/2+1/4 „ yoo y-j/2-3/4 ^ _ yoo yj/2-1 yoo yj/2-1 

X — Zj^-^qCjI i , (/ — Zjj-^gCj2 J- , 2: — Zj^-^qCj3 ± , u; — Zjj-^qCj4 ± , 

we arrive at the following asymptotic expansions: 

X = «Ti/^ + «C/3T5/' + C3iT^/^ + . . . , ^3 4^) 

y = a/4T-3/^ + 5aC/12Ti/^ + 7/4c3iT=^/^ + . . . , (3.42) 

z = v^/{A^/AX)T~^ -CV3/{3\^) -C3iV3/{a\^)T^/^ + (3.43) 

w = CT-^ - 4CV3 - 4Cc3i/aTi/2 + . . . . (3.44) 

The next step in our analysis is to check if for each j satisfying j/2 = g with g a positive 
eigenvalue, the corresponding eigenvector v of the transpose of the //Ci matrix is such 
that the compatibility conditions hold, namely, 

v'^ ■ Pj = 0, (3.45) 

where Pj are polynomials in Ci, . . . , Cj_i given by 

{lie, - (j72)X)c, = P,-. (3.46) 

Here the corresponding relation j/2 = 3/2, is valid only for j = 3 and the compatibility 
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condition indeed holds since 



(,/Ci-3/2X)c3 



V 



7 






1 


~4 




3 


3 


~16 


~4 






Aay/AX 




C 




a 


a 



a 



3 
2 







3 

'2 / 



C31 
7 



aVAX 



C31 



v 



4C 
C31 



(3.47) 



This means that a representation of the solution asymptotically with a Puiseux series 
as this is given by Eqs. fl3.4ip -( l3^4^ is valid. We thus conclude that as y — )■ F^, or 
equivalently as T — )■ 0, 



a — )■ 0, a' — )■ 00, 0' — )■ 00, P2 — ^ C)0. 
We consider also the case 7 = —1/2 for which we find 



(3.48) 



X = aT^/^ - aVlA/VScasT^/^ + . . . , 

y = a/4T-3/4 - 7ayiA/(4y3)c33T=^/^ + . 

z = y3/(4yiA)T-i + C33TI/2 ^ . . . , 

^ = ^ 2CyAA/V3c33T + . . . , 



(3.49) 
(3.50) 
(3.51) 
(3.52) 



where ( = —3kH'^/{2Aa'^). The compatibility condition for j = 3 is satisfied in the 
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following way: 



(,/Ci-3/2X)c3 



7 

4 



v 



3 


3 


~16 




73 








C 


-2<- 


2a 


a 



a 



V3AX 



7 VAXa 



3 
2 



3 

'2 / 



'4 V3 



-C33 



v 



C33 



2C^^C33 



73 



(3.53) 



We therefore see that for 7 = — 1/2 Eqs f l3.49p -( l332l) express the asymptotic behavior 
that corresponds to a curved brane around a collapse I singularity with the density of 
the fluid diverging, i.e as T — )■ 0, 



a — )■ 0, a' — )■ 00, 0' —7- 00, P2 — )■ 00. 

For 7 = —2 we find a different behavior 

X = aT^/^ - yiAa/73c33T^/^ + . . . , 

2/ = a/4T-3/^ - 7yiAa/(4y3)c33T3/4 + . 

^ = V3/(4yiA)T-i + C33TI/2 ^ . . . , 

w = CT - 4CVlA/y3c33T5/' + . . . . 



(3.54) 



(3.55) 
(3.56) 
(3.57) 
(3.58) 



We ought to check the compatibility condition for j = 3. We find that this is trivially 
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satisfied since 



(7/Ci-3/2X)c3 



_7 

~4 

3 

"16 



3 



a 



iaVAX aVAX 



3 
2 



V 



a 



a 




(3.59) 



The relations f l3.55p -f l3.58p are therefore valid representations of a general solution around 
the singularity at Y^. We can therefore conclude that as T — )■ 0, 



a — )■ 0, a' — )■ oo, — )■ oo, p2 — )■ 0. 



(3.60) 



A yet different behavior is met for 7 



-1: 



X = aT^/^ - yiAa/VScasT^/^ + . . . , 

y = «/4T-3/4 - 7yiAa/(4V3)c33T3/4 + . 

z = v^/(4yiA)T-i + ^ . . . , 

= ( + •••■ 



(3.61) 
(3.62) 
(3.63) 
(3.64) 
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The compatibility condition at j = 3 holds true since we find 



(,/Ci-3/2I)c3 



7 

1 



V 



3 


3 




~16 


~4 


2 






3 


Aay/AX 




~2 












3 

'2 / 



AXa \ 

C33 



7\/AXa 
C33 



-C33 



V 





V"/ 

(3.65) 



(3.66) 



so that the relations f l3.6ip - p.64p affirm that as T — )■ 0, 

a — )■ 0, a' — )■ 00, 0' — 00, p2 — )• C- 

3.4.2 The balance m-Bg 

We now move on to examine the next balance mBs = {{a, a/4, V^/i4:y/AX), 0), (1/4, -3/4, -1, -2)} 
which corresponds to a potentially general solution of a ffat brane. The A^- matrix for 
this balance is 

HiJCs = Df{a, a/4, V3/{aVAX), 0) - diag(l/4, -3/4, -1, -2), (3.67) 



where Df is the Jacobian matrix of the vector field f in Eq. (13. 5p . The eigenvalues of 
the ////C3 matrix are 

speciujlCs) = {-1, 0, 3/2, 1 - 7}. (3.68) 

Setting 7 = we get 

spec(m/C3) = {-1,0,3/2,1}. (3.69) 

We therefore have three non-negative /C-exponents which means that in this case the 
balance indeed corresponds to a general solution. The asymptotic expansions of the 
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variables in this case are 

X = aT^/^ + 2/3Aac2AT^^^-a^/AX/3c33T'^/^ + ..., 

y 



a/4T-3/4 + 5/6Aac24Ti/^ - 7/4a v^AA/ScgsT^/^ + . . . . 
z = y3/(4yiA)T-i - 2vW(3A)c24 + cssT^/' + • • • , 

W = C24T~^ + 



(3.70) 
(3.71) 
(3.72) 
(3.73) 



We have to check the compatibihty conditions for j = 2 and j = 3. Since we find 



(m/C3-X)c2 



and 



[m/C3-3/2X)c3 



5 
~4 


1 
i 


n 
U 


3 


1 


aV^AX 


^16 


4 


2 






-1 


Wax 


a VAX 











7 

~4 


1 





3 


3 


aV^AX 


~16 


~4 


2 


x/3 


Vs 


3 


4aA/AX 


ay/AX 


~2 
















{2/3)Aa 



\ / 2Aa \ 



hAa 
2y/A 



C24 



C24 



V C24 / 





(3.74) 



-(2/3)Aa 



1 
2 



V 



V3 

7 VAXa 
4 

C33 





C33 



-C33 



(3.75) 

the compatibihty conditions for j = 2 and j = 3 hold true. Eqs. fl3.70p -f lX75|) then imply 
that as T — 0, 

a — )■ 0, a' — )■ 00, 0' — c>o, p2 — ^ oo- (3.76) 
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In order to be able to compare and contrast our results found in this Subsection with 
the ones that are presented later in Subsection 3.7, we find it necessary to analyze here 
two more values of 7, namely 7 = —1/2 and 7 = —3/4. For 7 = —1/2 we find that the 
eigenvalues of the ////C3 matrix read 



speciniJC-i) = {-1,0,3/2,3/2}. 



For this value of 7 we find the asymptotic behavior 



X 

y 

z 
w 



a/4T-3/4 _ 7v/:4A/(4V3)ac33T='/^ + . 
C34T-V2 + . . . . 



Since 



(,,,/C3-3/2X)c3 



7 

~4 

3 

"16 



3 
"4 



a 



\ / VAXa \ 
7y/AXa 



V3 



4aVAX aVAX 



3 
'2 



V 











0/ 



4^3 



"C33 



V 



C33 



C34 



(3.77) 

(3.78) 
(3.79) 
(3.80) 
(3.81) 



/o^ 





(3.82) 

the compatibility condition for j = 3 is trivially satisfied, and as it follows from Eqs. 

dSZHD-dOT) as T ^ 0, 



a — )■ 0, a' — )■ 00, 0' — > c>o, p2 — )• 00. 



(3.83) 



Finally, for 7 = —3/4 the eigenvalues become 



speciniJC-i) = {-1,0,3/2,7/4}, 



(3.84) 
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while the asymptotic behavior now is 



X = aT^/^-VAX/V3ac63T^/^ + 16/33Aac74T^ + ..., (3.85) 

y = a/AT-^/^-7VAX/{AV3)ace3T^/^ + 32/33Aac74r + ... (3.86) 

z = ^/(4yiA)T-^ + cesT^/^ - 16y3A/(33yA)c74T3/^ + . . . , (3.87) 

w = C74T-^/^ + . . . . (3.88) 

We check the compatibihty conditions for j = 6 (note that here S = 4 and we see that for 

the eigenvalue g = 3/2 the corresponding arbitrary constant appears when j = gS = 6) 
and for j — 7 are again trivially satisfied since we find that: 



(m/C3-3/2X)c6 





1 


4 




3 


3 


"16 


~4 


V3 






ay/AX 



and 



(777/C3-7/4X)c7 



3 

"16 



-1 



4q;\/A\ aVAX 








3 
2 



a 





v/3]4A 



7 
'i 







1/3Aq; 



MA 



4\/3 



aC63 



-ttC63 



C63 



1/4 / 






















l/3Aa 



16 , 
33 ^ 



32 , 

— Aac7A 

33 ^ 



(3.89) 



16 



-^/3A/Xcu 



33 



J V 



/ V ° / 

(3.90) 
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Hence Eqs. flSlSSjl - flSlSH]) show that as T ^ 0, 

a — )■ 0, a' — )■ cxD, 0' — )■ oo, p2 — )■ oo. 



(3.91) 



3.4.3 The balance jB2 



Here we consider the balance 1B2 that corresponds to a fiat brane. The eigenvalues of 
the /C-matrix for this balance are 



spec(7/C2) = {-1, 0, 3/2, -1-7-4. 
For s = 1/2 and 7 = —5/2 we get 

spec(7/C2) = {-1,0,3/2,1}. 
We then find the following expansions 



X 

y 

z 
w 



a 



«/4T-3/4 _ 7/4V'Ia73«C33T3/^ + 

C24T='/2 + . . . 



We have to check the compatibility conditions for j = 2 and j = 3. Since 

1 \ 



4 



[//C2-X)C2 



3 


1 


^16 


^4 

















VSAX 



-1 



/ A 



/ 



y C24 y 



(3.92) 



(3.93) 



(3.94) 
(3.95) 
(3.96) 
(3.97) 



(3.98) 
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and 



(,/C2-3/2X)c3 



7 

1 



V 



3 


3 




^16 




2 




Vs 


3 


4a a/A A 




~2 












-C33 







1 

2 / 



4 73 



C33 



V 











(3.99) 

the compatibility conditions are indeed satisfied and from Eqs. fl3.94p - (l3.97p we see tliat 
as T ^ 0, 

a 0, a' 00, </)'-)■ 00, p2 0. (3.100) 



3.4.4 The balance jjjBi 

The balance mBi that corresponds to a potentially general solution of a flat brane is 
valid only for 7 = 1 and describes the asymptotics around a collapse type I singularity. 
Here the eigenvalues of the matrix 777/Ci are 

spec(m/Ci) = {-1, 0, 0, 3/2}. (3.101) 

Since we find three non-negative /C-exponents for this balance we conclude that it indeed 
corresponds to a general solution. 

Here we find the following expressions 

X = aT^/^ - «/(4c)c33T^/^ + . . . (3.102) 

y = a/4T-3/^ - 7a/(16c)c33T=^/^ + . . . (3.103) 

z = cT-i + C33TI/' + . . . (3.104) 

w = (3/(32A)-AcV2)T-2 + (2/c)(3/(32A)-AcV2)c33T-i/2 + ... (3.105) 
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and the compatibility condition for j = 3 is valid since 
/ 7 



(m/Ci-3/2X)c3 



4 
3 

"l6 

c 
a 



1 



' 'Ac''' ^ 



3 

4 



-4- 



\ 16ayl a 



a 



3 4Ac2 
+ 



-2A\ca -2Aa 

-2 

2 



4aA 







3 
2 



7a 
16c 



C33 



C33 



v 



-CC33 

c 



(3.106) 

where C = 3/(32^) - Xc'^/2. We thus conclude that as T ^ 0, Eqs. fl3302D -f l3J05D 
describe the asymptotics around a collapse I singularity with the density of the fluid 
diverging there, i.e., 



a — )■ 0, a' — )■ 00, — )■ 00, p2 — )■ 00. 



(3.107) 



3.4.5 The balance 



Finally, the balance 7/7-62 corresponds to a potentially general solution of a flat brane 
and has the following /C-exponents 



spec(777/C2) = {-1,0, -2(p - 1), 1 - 4p}. 



(3.108) 



The last two /C-exponents are positive when 7<— lor7>l. We consider here the 
case 7 > 1 because as we show below it leads to the emergence of a collapse I singularity 
while the case 7 < — 1 implies the existence of a big rip singularity and it is considered 
later in Subsection 3.6. We set 7 = 2. Then p = 1/6 for which the balance and the 
/C-exponents read 



and 



111B2 = {{a, a/6, 0, 1/(24A)), (1/6, -5/6, -1, -2)} 
spec(777/C2) = {-1,0,5/3,1/3}. 



(3.109) 
(3.110) 
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Since we have three non-negative /C-exponents we see that this balance indeed corre- 
sponds to a general solution. The variables in this case expand as follows, 



X 

y 

z 
w 



— a 



^ 3/5q,^a4t5/6 - 27/14q;>1'A24t3/2 - 2q;Ac54T^^/^ + . . . (3.111) 
^ i/2q;^a4T-^/6 - 81/28q;^'A'4T^/' - U/SaAc^^T^^^ + ■ ■ ■ (3.112) 

C13T-2/3 - 12/5AAc?3 + Sge/SSA^A^c^gT^/^ + . . . (3.113) 
l/(24A)T-2 - 3/10Ac?3T-^/3 ^ 747/350AX^4^T-^/^ + 05^^'^^^ + ... (3.114) 



We ought to check the compatibility conditions for j = I and j = 5. We find 

[ -1/2 1 

-5/36 1/2 -10/3aA 

ci3 

\ l/(12aA) -l/{2aA) -1/3 J \ ^ / 



(777/C2 - l/3X)ci 



(3.115) 

so that the compatibility condition for j = 1 is satisfied. Also for j = 5 we have 



(m/C2 - 5/3X)c5 



/ -11/6 


1 





\ 


-5/36 


-5/6 





-10/3q;^ 








-4/3 





^ l/{12aA) 


-l/(2a^) 





-5/3 ; 


( ° 1 










-4/3C53 









/ -2aAcu \ 
— 11/3q;Ac54 

C53 
C54 



(3.116) 



V / 

where C53 = 396/35 A^A^cf 3. The corresponding eigenvector v is such that 

v'^ = (l/(12/la),-l/(2/la),0,l). 
The compatibility condition, 

• Pj = 0, 



(3.117) 
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for j = 5 therefore holds true and Eqs. fl3.llip - fl3.114p represent the asymptotics around 
a collapse I singularity with the density of the fluid being divergent, i.e. as T — )• 

a — )■ 0, a' — )■ oo, 0' — c>o, p2 oo. (3.118) 
3.5 Collapse type II singularity 

We shall analyze in this Section the asymptotics represented by the balances uB^ and 
11184^ which suggest the emergence of a collapse type II singularity. 

We start with the balance 11B3 that corresponds to a curved brane. The /C-matrix 
for this balance is given by 

J7/C3 = L'f(°na,a,0,0) -diag(l,0,r,-2), (3.119) 

where Di^^^ is the Jacobian matrix of the dominant part f*-"-* in Eq. f l3.2ip . The eigen- 
values of ///C3 are 

spec(,,/C3) = {-1, 0, -2 - 47, -4 - r}. (3.120) 
For 7 = — 1 and r = 1, this balance becomes 

7/^3 = {(a, a, 0,0), (1,0, 1,-2)}, (3.121) 

and the eigenvalues of the A^-matrix are 

spec(,,/C3) = {-l,0,2,-5}. (3.122) 

We may set the arbitrary constant appearing for j = —5 equal to zero and find a 
particular solution around a collapse II singularity. The asymptotic expansions of the 
variables now are: 

X = aT + Aa/9c24T^ + ... (3.123) 

y = a + Aa/3c24T^ + ... (3.124) 

z = + ... (3.125) 

w = C24 + ... . (3.126) 
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(///Ca - 2X)C2 



The compatibility condition for j = 2 is satisfied since 

/ -3 1 \ / (aA/9)c24 \ 
0-2 (2/3)aA (aA/3)c24 
0-30 
V / \ C24 y 
It follows then from Eqs. (I3323|) - (13J26D that as T ^ 0, 

a — )■ 0, a' — )■ a, 0' — 0, p2 — ^ C24. 



/o\ 






(3.127) 



(3.128) 



Next we turn to the balance mB^. This balance corresponds to a potentially general 
solution of a curved brane. The eigenvalues of its /C-matrix are 



spec(m/C4) = {-3, -2 - 47, -1, 0}. 



(3.129) 



For 7 < —1/2 we have one positive eigenvalue (the case 7 > —1/2 for which we have 
three negative eigenvalues is considered later in Subsection 3.7). We may set c_3 = 
from the beginning and find an asymptotic expansion of a particular solution around a 
collapse II singularity. We choose 7 = —3/4. Then 



spec(m/C4) = {-3, 1, -1, 0}, 
and we find the following asymptotic forms 

X = a 

y = a + Aa/3ci4T + . . . 
z = + ... 

W = Ci4T~^ + . . . . 

The compatibility condition for j = 1 is trivially satisfied since 

/ -2 1 \ f (aA/6)ci4 \ 



(3.130) 

(3.131) 
(3.132) 
(3.133) 
(3.134) 



[1 1 1 '^4 



/C4 — X)ci 



0-10 aA/3 
0-40 




{aA/3)cii 


Ci4 



/o\ 






(3.135) 
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The forms in Eqs. f l3.13ip - fl3.134p then imply that as T — )■ 0, 

a 0, a' a, 0' 0, p2 oo. (3.136) 

Thus for a curved brane we found collapse II singularities with 0' — )■ and p2 — ?■ oo, ps, 
Ps ^ 0. This means that asymptotically the leak of energy from the brane and into the 
bulk is controlled solely by the fluid. 

3.6 Big rip singularity 

In this Section we examine the balances uBi and 111^2 which correspond to a flat brane 
and describe the asymptotics around a big rip singularity when we choose a value of 7 
such that 7 < — 1. 

The first of these balances has the following /C-exponents 

spec(,,/Ci) = {-1, 0, 0, ~2{p - 1)}. (3.137) 

Choosing 7 = —2, we have that 

iiB, = {{a, -a/2, c, 3/(8A)), (-1/2, -3/2, 2, -2)} (3.138) 

and 

= {-1,0,0,3}. (3.139) 

This leads to the following solution, 

X = aT-^/^-a/{4c)c3^T'/^ + ... (3.140) 

y = -a/2T-3/2 - 5a/(8c)c33T=^/^ + . . . (3.141) 

z = cT^ + C33T5 + . . . (3.142) 

w = 3/(8A)T-2 - 3/(8Ac)c33T + . . . (3.143) 
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We validate the compatibility condition for j = 3, since 



(///Ci - 3X)c3 



/ -5/2 


1 





\ 


/ 


—Q:/(Ac)C'i'^ ^ 








3/4 


-3/2 





2aA 




-5a/(8c)c33 









-2c/a 


-Ac/a 


-3 
















^ 3/(4aA) 


3/{2aA) 







v 


-3/(8Ac)c33 yl 






















(3.144) 


a big rip sinj 


jularity: 
















a — oo, 


a' — )■ — oo 






P2 


— )■ CXD. 




(3.145) 



We now examine the second balance 11182- The /C-exponents for this balance are 
given by Eq. f l3.108p . We are interested in the case 7 < — 1 which we have not examined 
yet and for which we get two positive /C-exponents. As we show below this case leads to 
the emergence of a big rip singularity. We set for concreteness 7 = —3/2. Then 



jnB2 = {{a, -a, 0, 3/(2A)), (-1, -2, -1, -2)} 



(3.146) 



and 



////C2 = {-1,0,4,5}. 



(3.147) 



In this case we have three non-negative /C-exponents so that this balance indeed corre- 
sponds to a general solution. The variables in this case expand as follows, 



X = aJ-^ + Aa/3cuT^ + . . . 

y = —aT^^ + Aac^iT"^ + . . . 

Z = C53T^ + . . . 

w = 3/(2A)T-2 + C44T2 + . . . 



(3.148) 
(3.149) 
(3.150) 
(3.151) 
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We ought to check the compatibihty conditions for j = A and j = 5. We find 

/ -3 1 ^ \ I ^a/3c44 \ / \ 

2 -2 (4/3)aA 

0-10 

\ 3/(aA) 3/(aA) -4 



[111IC2 - 4X)C4 



Aac44 


C44 








(3.152) 



and 



(///^2 - 5X)C5 



-4 











\ / \ 



C53 
V / 



/o\ 






(3.153) 



1 

2 -3 (4/3)aA 



\^ 3/{aA) 3/{aA) -5 

so that Eqs. fl3.148p - fl3.151l) represent the asymptotics around a big rip singularity i.e. 
as T ^ 

a ^ 00, a' — > —00, 0' — > 0, p2 — c>o. (3.154) 

Hence in both of the cases studied in this Section we found that for a fiat brane, a big 
rip singularity develops at a finite distance. We note here that the exchange of energy 
from the brane into the bulk is totally monitored by the fiuid. 



3.7 Behavior at infinity 

In this Section we consider the balances 1B3, mB^, jB^, 11B2 and inB2 that offer the pos- 
sibility of escaping the finite-distance singularities met before and describe the behavior 
of our model at infinity. 

We begin with the analysis of the balance iB-^. The eigenvalues of its /C- matrix read 

spec(,/C3) = {-1,-3,0,0}, (3.155) 

hence we may expand (x, w) in descending powers in order to meet the arbitrary 
constants appearing at j = —1 and j = —3. We choose 7 = —1/2. The balance jB^ then 



30 



corresponds to a general solution of a flat or curved brane. In particular we find 



X = aT + c_ii + . . . 
y = a + ... 



(3.156) 
(3.157) 
(3.158) 



w = CT"^ - 2C/ac„iiT-^ + 3C/aV uT"^ - 4C/aV nT"' + . . . . (3.159) 
The compatibility conditions for j = —1 are satisfied since 



(7/C3+X)c_i 



For j = — 3 we find 
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C-11 
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V 2C/« 
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-2C/«c 


























(3.160) 



(j/C3+3X)c_ 



( ^ 


1 








\ 




\ 
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C-33 









V 2C/« 


-2C/a 
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\ 




J 



P. 



3- 



Since a corresponding eigenvector here is v'^ = (0, 0, 1, 0) and 

v'^ ■ P_3 = 



(3.161) 



(3.162) 



the compatibility condition for j = —3 is also satisfied and Eqs. (I3.156I) - (I3.159I) then 
imply that as T — )■ oo 



a — )■ oo, a' — )■ a, </>'—)■ 0, p2 0. 



(3.163) 



We now examine the balance hiISa = {{a, a, 0, 0), (1, 0,-1, —2)} that corresponds to 
a general solution of a curved brane. The /C-exponents are given by Eq. (13.129^ . For 
7 > —1/2 we have three negative /C-exponents. We choose 7 = 0. Then 



spec(////C4) = {-3, -2, -1, 0}. 



(3.164) 
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We find here that 



X = aT + c_ii - (1/3)q;Ac_24T-^ - (5/9)Ac_iiC_24T-^ + 

y = o; + (l/3)Q;Ac_24T-' + (10/9)Ac_nC_24T-^ + --- 

z = c_33T-^ + ... 

w = c_24T-^ + (4/a)c_nC_24T-5 + . . . . 



(m^4 +2:)c_i 



while for j — —2 



{iii^i + 2X)c_2 
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-(2/3)q;^ 
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1 1 







v 



2 -(2/3)aA 
0-1 




( -(l/3)aAc_24 \ 
(l/3)aylc_24 


V ^-24 J 



/o\ 






For j = — 3 we have that 



(7/7^4 + 3X)c_3 



/ 2 1 \ 



3 -(2/3)q;^ 

1 

/ \ 

2/3Ac_iic_24 


y 4/q;c_iiC_24 y 



( -(5/9)Ac_nC_24 \ 

(10/9)>1C_11C_24 
C-33 

y (4/q;)c_iiC_24 y 



(3.165) 
(3.166) 
(3.167) 
(3.168) 



We check the compatibility conditions for j = —1, j = —2 and j = —3. For j = —1 we 
find 

/o\ 



(3.169) 



(3.170) 



(3.171) 
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and an eigenvector v is such that = (0, 0, 1, 0). The compatibihty condition, 



V'^ ■ P-3 = 0, 



(3.172) 



therefore holds true. Eqs. fl3.165p -( T3.168p then imply that as T — )• oo 

a — )■ oo, a' — )■ a, 0' — 0, p2 — )■ 0. 



(3.173) 



On the other hand, the /C-exponents for the balance 1B4, that also corresponds to a 
curved brane, are 

spec(7/C4) = {-3, -s - 4(7 + 1), -1, 0}. (3.174) 
For s = —3 and 7 = 1/4 we find the following behavior 



X 

y 

z 
w 



aT + c_ii - (l/6)aAc_24T'2 + 
a + (l/3)aAc_24T"^ + ... 
C-33T-^ + . . . 



-33 

-24 ■ 



c_,4T-5 - (5/a)c„nC_24T-'' + 



(3.175) 
(3.176) 
(3.177) 
(3.178) 



The compatibility conditions for j = —1 and j = —2 are trivially satisfied since 



(,/C4+X)c_i 








1 








\ 




C-11 


\ 







1 

























-2 














v 











-1 


/ 


v 





/ 



and 



[//C4 + 2X)c_2 



110 



2 
0-10 
\ J 



( \ 










/o\ 






(3.179) 



(3.180) 
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For j = —3 we have that 



(7/C4 + 3X)c. 



2 10 



/ -(l/6)Aac_24 \ 



(l/3)Aac_24 
C-33 

\^ -(5/a)c_iic_24 / 
\ 



(3.181) 



3 


\ 1 y 

/ 

Aac-2i 


y -(5/a)c_iic_24 / 

and an eigenvector v is such that = (0, 0, 1, 0) so that the compatibihty condition for 
J = -3, 

yT.p_3 = o, (3.182) 
is also satisfied. Therefore Eqs. fl3.175p -( [3.178p then imply that as T — )• 00 



a — )■ 00, a' — )■ a, 0' — )■ 0, p2 — ^ 0. 



(3.183) 



Let us now examine the balance 11B2 that corresponds to a flat brane. For 7 = —3/4 
(hence p = 2) and r = — 5 this balance reads 11B2 = {(a, 2a, 0, 6/ A), (2, 1, —5, —2)}. and 
the eigenvalues of the 7-7-/C2 matrix are 



spec(///C2) = {-3, -2, -1, 0}. 



(3.184) 



We find here the following asymptotic behavior 



X 

y 

z 
w 



aT^ — Aa/Gc^i^T — y4a;/6c_24 + A^a/SGc^^^ + . . . 
2aT — Aa/6c^i4 + . . . 



(3.185) 
(3.186) 

-33^ (3.187) 
+ c_i4T-3 + c_24T-^ + (-AV36c^4 + A/3C_14C_24)T-^ + .(3.188) 
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The compatibility conditions ior j — —1, j — —2 and j — —3 are satisfied. Particularly 
for j' = — 1 we find 

/ -1 1 \ / -Aa/6c_i4 \ / \ 

2 Aa/3 

0-2 

\ 12/{Aa) -6/{Aa) 1 J 



(,,/C2+X)c_i 



—Aa/6c-i4 


V ^-14 / 








For j — —2 we find that 
/ 



(/J/C2+2X)C_2 



1 \ / - 

2 1 Aa/3 

0-10 

\^ 12/(Aq;) -6/{Aa) 2 J 

( \ 

A^ajl^c^^^^ 




(3.189) 



2 . \ 



\ / -^a/6c_24 + ^V36di4 



C-24 



V 



(3.190) 



and an eigenvector v is such that 

= (12/(Aq;), -6/(Aq;), 0, 1) 



(3.191) 



and hence we find 

■ P_2 = 0, (3.192) 

which means that the compatibihty condition is satisfied also for j = —2. Finally, for 
J = — 3 we see that 

/ 1 1 \ / \ 

2 2 Aa/3 

C_33 

\l2/{Aa) -6/{Aa) 3/ \ -^736^ 14 + A/3c_i4C_24 y/ 



(,,/C2+3X)c_3 
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/ \ 

Aa/3(-AV36c3 ^4 + A/3c_i4C_24) 


\ 3(-AV36c3,4 + A/3c_i4C_24) ) 
while an eigenvector v is such that 

t;T = (0,0,1,0), 



(3.193) 



and hence 



V ' ■ P-3 = 0, 



(3.194) 



(3.195) 



so that the compatibihty condition for j = —3 is satisfied. We thus see from Eqs. 
(13:T85D -( 13:T88D that as T ^ cx) 



a — )■ oo, a' — i- oo, 0' — 0, p2 — )■ 0. 



(3.196) 



We now move on to the balance 111^2 which also corresponds to a fiat brane. It follows 
from fl3.108p that for— 1<7<— 1/2 this balance leads to three negative /C-exponents. 
We take 7 = —3/4 and then this balance becomes 



inB2 = {{a, 2a, 0, 6/ A), (2, 1, -1, -2)}, 



and the eigenvalues of the ////C2 matrix are 



spec(m/C2) = {-7,-2,-1,0}. 



In this case we find that 



X = aT^ + c_iiT + c_2i + 
y = 2aT + c_ii + ... 



w = 6/AT-'-6/(Aa)c_iiT-^ + 6/(Aa')(ciii -ac_2i)T-4 + 

+ C_34T-^ + C_44T-6 + C_54T-' + C_64T-^ + 0.^4^'^ + • • • : 



(3.197) 

(3.198) 

(3.199) 
(3.200) 
(3.201) 

(3.202) 
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where Cj4 with j — —3, . . . , —7, are polynomials in a, C-n and c_2i. The compatibility 
condition at j = — 1 is trivially satisfied since we have that 



(m^2 +2^)c-i 



/ -1 
2 




1 






\ / 

Aa/3 
-6 



\^ 12/{Aa) -6/{Aa) 1 ^ 



C-ll 
C-ll 



\ -6/(Aq;)c_ii J 



For j — —2 we find that 
/ 



(m/C2+2X)c_2 



2 




1 \ / 

1 Aa/3 

-5 



12/{Aa) -6/{Aa) 2/ 



C-21 






/o\ 




(3.203) 
\ 



\^ 6/(Aa2^(dii -q;c_2i) / 

(3.204) 



2/acl 



11 



(3.205) 



V 12/ {Aa')cU, I 
and t;^ = {12 /{Aa), —Q/{Aa), 0, 1) so that the compatibility condition 

• P_2 = 

is satisfied. Finally, for j — —7 we have that 

f 5 1 Oo\/o\ 

2 6 ^a/S 

C_73 



(m/C2+7X)c_7 



/ \ 

ylQ;/3c_74 



12/(Aa) -6/{Aa) 7 J \ 0.74 / 







V 



7C_74 J 

(3.206) 

and the corresponding eigenvector here is such that v'^ — (0, 0, 1, 0) which implies that 
the compatibility condition 



- P. 
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holds true. From Eqs. fl3199|) - (l3:202|) it follows that as T ^ cx) 

a^oo, a' ^ oc, 0' ^ 0, p2 0. (3.207) 

In our previous work in [2] , the bulk was filled entirely by the fiuid and we found that 
for— 1<7<— 1/2 and for a fiat brane the avoidance of finite-distance singularities was 
the only possible asymptotic behavior. This fact suggested that a self-tuning mechanism 
could be build within the framework of such model, a property that we would anticipate 
to hold also in the more complicated case studied in this paper. However, by the analysis 
we did so far, we see instead that although the balance 7-//-B2 implies a behavior that 
is singular only at infinite distance, the balance mB^ that is also valid for this range 
of 7 implies a singular (collapse type I) behavior at finite distance (as this follows from 
Eqs. f l3.85l) - p.88p in Subsection 3.4.2.). Both of these balances represent behaviors of 
the general solution. In particular, the balance 111^2 describes the asymptotic behavior 
of the general solution in a neighborhood of infinity, while, the balance 7/7-83 describes 
the asymptotic behavior of the general solution around a finite-distance singularity. We 
thus conclude that in this case the avoidance of finite-distance singularities for fiat brane 
becomes impossible. Our results suggest that the singular behavior encountered here 
is driven by the presence of the scalar field which is left to act independently from the 
fluid. What would happen, instead, if it interacted with the fiuid? In the next Section, we 
show that by choosing the interaction parameters in an adequate way, we may resolve this 
unwanted situation and recover the possibility of avoiding finite-distance singularities. 

4 Interacting mixture in the bulk 

In this Section, we study the possible behaviors that arise when the two bulk components 
interact with each other. We begin by searching to find what are the forms of the balances 
that are possible in this more intricate case. In order to simplify our calculations, that 
are much more complicated than in the case studied previously, we may set to zero either 
one of the two parameters a and u that define the interaction and let the remaining one 
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vary arbitrarily. If we choose both parameters nonzero we are led to balances {a, p} with 
the exponents in the vector p being irrational which leads to the existence of logarithms 
in the series expansions of the variables. In the next paragraphs we show that the choice 
that gives the desired result, meaning the avoidance of singularities, is a = oO 

We start the analysis by putting a = in the system f l2.11l) - fl2.14p and letting v be 
arbitrary but nonzero. We consider all terms dominant and by substituting Eq. fl3.10p 
there, we find the following four balances 




2a / 3(47 -4- z/) 3z/ 



8 + z/' V AA(7 - 1)(8 + z/)2' 2A(7 - 1)(8 + vf 




3x = 

I \ s -|- Z/ ' 

^62 = { (a,ap,0,^^ 1,-1,-2)|, p= 2(^^\i^ 

= {(a, a, 0,0), (1,0, -1,-2)}, 
= {(a, 0,0,0), (0,-1, -1,-2)}. 

The balance is valid for u 7^ —6, —8,47 — 4 and 7 7^ 1 and because of the square 
root we may have either > 47 — 4 and 7 < 1, or, z/ < 47 — 4 and 7 > 1. Substitution of 
these balances in the constraint equation ( 12.15^ shows that ^i^i, and 1,84 correspond 
to a flat brane whereas corresponds to a curved brane with a satisfying a^ = kH^. 



"^The case v — and a arbitrary does not lead to avoidance of singularities but instead it brings back 
the same problem we faced in Subsection 3.7. 
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We now calculate the Jacobian of the vector field (12.161) (with cr = 0) and find 



/ 1 \ 



--A{l + 2-f)w - XAz^ -2XAzx --A{l + 2-f)x 

3 3 

-f4+-)- 

2/ x2 V 2/ X V 2/ a; 

\ 4 7 + 1—-- ^ -4(7 + 1)- + - uX— -4(7 + 1)- 

\ 1 x'^ X 2 X X X ' 



(4.5) 



The balance yB^ is discarded because it does not give the —1 /C-exponent, but it has 
instead 

spec(,/C4) = {0,1,1,2}. (4.6) 
The balance j^Sa, on the other hand, has 

spec(,/C3) = {-1, 0, -2(1 + 27), -3 - Z//2} , (4.7) 

which implies that for a curved brane we may avoid the finite-distance singularities for 
7 > — 1/2 and v > —6. We may also include in this range of 7 the value —1/2 since then 
the following balance arises: 

,i3_i/2 = {(a, a, 0, 6), (1, 0, -1 - 2)}. (4.8) 

The balance i,B^i/2 corresponds to a general solution of a fiat or curved brane with 
6 = 3/{2A){l - knya^), (5 ^ and gives 

spec(,/C_i/2) = {-1, 0, 0, -3 - Z//2}, (4.9) 

so that for u > —6 finite-distance singularities can be avoided. 
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The /C-exponents for the balance i^Bi are 



. ,^ s f 2(6 + u) 2(4 - 47 + z/) 1 
.pecW^j-l.O.^. > y (4.10) 

while for the balance we find, 

spec(./C2) = 0, ^^^^ |. (4.11) 

The balance ,^^2 implies that we may escape finite-distance singularities for — 1 < 7 < 
— 1/2 and u > —4 + 47 since then the last two eigenvalues of u}(^2 become negative. We 
note though that for these ranges of 7 and u the balance i,Bi is also valid and it is such 
that at least the last eigenvalue of j//Ci becomes then positive. We may exploit this fact 
by keeping 7 in the interval (—1, —1/2) and restricting u to fall in (—4 + 47, ~6). For this 
choice of parameters, the last two eigenvalues of ^Ki have opposite signs which means 
that the solution described by is neither valid around a finite-distance singularity, 
nor, around a neighborhood of infinity but rather in the limited area of an annulus 
failing thus to provide us with any substantial information about the asymptotics of the 
dynamical system fl2.1ip -( l2?l^ . We are then left with the unique possibility described by 
the balance iyB2 which is the asymptotic expansion of the general solution of the system 
in a neighborhood of infinity. Consequently, avoidance of finite-distance singularities is 
feasible for -1 < 7 < -1/2 and -4 + 47 < z/ < -6. 



5 Conclusions 

We have studied a model consisting of a three-brane embedded in a five-dimensional bulk 
filled with a scalar field and an analog to a perfect fluid possessing a general equation 
of state P2 = 7P2; characterized by the constant parameter 7. The two bulk matter 
components may act independently, or, they may interact with each other by exchanging 
energy in a way that the total energy is conserved. 

We have started off by analyzing the evolution of our model in the case that the two 
bulk matter components do not interact with each other (the two interaction parameters 
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are set to zero in this case) . We have shown that a quite general feature of the asymptotic 
behavior of such model is the emergence of a finite-distance singularity that is of the 
collapse type I, II or big rip class. The singularities accommodated in the first two 
classes share common characteristics such as the vanishing of the warp factor. However, 
the derivative of the warp factor behaves differently in each case: it is divergent in the 
collapse type I class whereas it remains finite in the collapse type II class. The collapse 
type I singularity may arise for all values of 7, whereas, the collapse type II class arises 
only when 7 < —1/2 as this is illustrated by the balance 111B4 but also by the balances 
11B3, jBs and 7B4 for an adequate choice of their parameters. All of these balances lead 
to particular solutions for this range of 7. The third class, on the other hand, of big rip 
singularities, arises always with 7 < — 1 and it is characterized by the divergence of the 
warp factor, its derivative and density of the fluid while the energy density of the scalar 
field now tends to zero in the neighborhood of the singularity. 

We completed our analysis for non-interacting bulk matter by addressing the impor- 
tant issue of whether it is possible to avoid finite-distance singularities. We found that 
this is true only for a curved brane and for 7 > —1/2. This is demonstrated by the bal- 
ances 111B4, jjBs and /B4 that give negative /C-exponents for 7 > —1/2. For a flat brane 
on the other hand, the avoidance of singularities is not possible, as this was discussed in 
detail in the last paragraph of Subsection 3.7. The main reason for the failure to escape 
finite-distance singularities in the case of flat brane is that for all values of 7 there al- 
ways exists a balance that corresponds to a general solution and describes its behavior 
around finite-distance singularities due to the presence of the scalar field component in 
the bulk. It is thus impossible to find ranges of 7 that they are not characterized by 
singular behavior. 

For illustration, we present a summary of our results for this first analysis in the 
Table 1 below, using the notation for the various singularities introduced in Section 2 
and the symbol * to denote a balance that corresponds to a particular solution. For 
each entry in the table we have taken into account all the corresponding examples in our 
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analysis to deduce the form of each balance and /C-exponents. Note that we used the 
numerical examples only as representatives of the corresponding asymptotic behaviors for 
the different regions of the parameter 7. We also give the ranges of r and s, entering in the 
definitions of the balances around the type I and II singularities defined in Subsections 
3.1 and 3.2, that lead to the most general behavior possible for each balance. 

We continued our analysis to include also cases with interacting bulk matter, moti- 
vated by the fact that it is impossible to find in our flat brane model ranges of 7 that 
lead to avoidance of singularities within finite distance in the case of non-interacting 
bulk matter. We studied the case of interaction a — 0, u arbitrary. We have shown 
that this choice leads to the avoidance of finite-distance singularities for— 1<7<— 1/2 
and —4 + < u < —6. For a curved brane, avoidance of finite-distance singularities is 
allowed for 7 > —1/2, u > —6 and a = 0. These results enforce our previous conclusion 
about the possibility of such solutions and show that they are also possible in the more 
involved case considered in this paper, thus proving that the self-tuning mechanism may 
be sustained under field interactions in the bulk. 

We illustrate the results we found for the case of non-interacting bulk components as 
well as for the case with interaction a — and ly arbitrary in the following tables: 
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equation of state 


flat brane 


curved brane 




type 


balance 


type 


balance 


7 > 1 


singular type I 


/jSi, 111B2, 
11B2 (-l<r< -2/(7 + 1)) 


regular 


at 00: 111B4, 
1B4. (-4(7+1) <s< -2), 
11B3 (-4 < r < -1) 


7=1 


singular type I 


///Si 


regular 


at 00: 111B4,, 
1B4 (-8 < s < -2), 
iiBs (-4 < r < -1) 


-l/2<7< 1 


singular type I 


iBi, iiiBs, 
1B2 (-2<s<-(l+7)) 


regular 


at 00: 111B4, 
iBi (-4(7+1) <s< -2), 
11B3 (-4 < r < -1) 


7 = -1/2 


regular 
singular type I 


at 00 1B3 
1B2 {s < -1/2), 111B3 


regular 
singular type I 


at 00: /S3, JB4 (s > —2) 
iB, 


-l<7<-l/2 


regular 
singular type I 


at 00: 111B2, 
11B2 (-2/(7+1) <r<-l) 
iBi, 111B3, 
1B2 (-2<s<-(l + 7)) 


regular 
singular type II 


iiiBl iBl {s < -2), 
iiB; (-4 < r < -1) 
jBl (-2<,s<-4(l+7)), 
iiiBl iBl, iiBl (r > -1) 


7 = -l 


singular type I 


iBi, 111B3, 
1B2 (-2 < s < 0) 


regular 
singular type II 


at 00: iBl {s < -2), iiiBl, 

iiBl (-4 < r < -1) 
iBl (-2<s< -4(1 + 7)), 
niBl, iBl iiBt (r > -1) 


7 < -1 


singular big rip 
singular type I 


iiBi, 111B2, 
11B2 (-l<r< -2/(7+1)) 
iBi, iiiBs, 
1B2 (-2<s<-(l+7)) 


regular 
singular type II 


at 00: iBl (s < -2), iiiBl, 

iiBl (-4 < r < -1) 
iBl (-2<s< -4(1 + 7)), 
niBl, iBl nB*3 (r > -1) 



Table 1: Summary of our results for the case of non- interacting bulk components. 
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equation of state 


flat brane 


curved brane 


P2 = 7P2 


type 


balance 


type 


balance 


7 > 1 


regular 
singular big rip 
singular type I 


lyBi: — 8 < < — 6 at 00 
uBi: u <-8 
yB2: z/ < -4 + 47 


regular 


,yB^: V > —6 at 00 


7 = 1 


singular type I 


uB2: < 


regular 


1^63: v > —6 at 00 


-1/2 < 7 < 1 


singular type I 


uB2: V < -4 + 47, 
^Br. z/ > -4 + 47 


regular 


lyB^: v > —6 at 00 


7 = -1/2 


regular 
singular type II 
singular type I 


uB_i/2- > —6 at 00 
„B_i/2- v <-6 
^Bi: -6 < 1/ < 


regular 
singular type II 


vB_i/2' ^ > ~6 at 00 
vB-il2- z/ < -6 


-1< 7 < -1/2 


regular 
singular type I 


i,B2: -4 + 47 < z/ < -6, 
u > —6 at 00 
^Bi: iy> -6 


singular type II 


^Hg: z/ < -6 


7 = -l 


singular type I 


^Bi: -6 < 1/ < 


singular type II 


^Bg: z/ < -6 


7 < -1 


singular big rip 
singular type I 


z/ > -4 + 47 
^Bi: u> -6 


singular type II 


^Hg: z/ < -6 



Table 2: Summary of our results for the case of interacting bulk components. 
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